Presented as a Keynote Lecture at the 15th U.S. National Congress of 
Theoretical and Applied Mechanics in Boulder, Colorado, on June 26, 
2006, and discussed in the Symposium on "Revisiting the Collapse of the 
World Trade Towers " 



Mechanics of Progressive Collapse: 

Learning from World Trade Center 
and Building Demolitions 



Zdenek P. Bazant and Mathieu Verdure 



Report No. 06-06/C605t 



Department of Civil and Environmental Engineering 

Northwestern University 

Evanston, Illinois 60208, USA 



June 23, 2006 



Submitted on June 23, 2006 

Mechanics of Progressive Collapse: 

Learning from World Trade Center and Building Demolitions 

Zdenek P. Bazant 1 , F. ASCE, and Mathieu Verdure 2 

Abstract: Progressive collapse is a failure mode of great concern for tall buildings, and is also 
typical of building demolitions. The most infamous paradigm is the collapse of World Trade Center 
towers. After reviewing the mechanics of their collapse, the motion during the crushing of one 
floor (or group of floors) and its energetics are analyzed, and a one-dimensional continuum model 
of progressive collapse is developed. Rather than using classical homogenization, it is found more 
effective to characterize the continuum by en energetically equivalent snap-through. The collapse, 
in which two phases — crush-down followed by crush-up — must be distinguished, is described in each 
phase by a nonlinear second-order differential equation for the propagation of the crushing front of 
a compacted block of accreting mass. Expressions for consistent energy potentials are formulated 
and an exact analytical solution of a special case is given. It is shown that progressive collapse will 
be triggered if the total (internal) energy loss during the crushing of one story (equal to the energy 
dissipated by the complete crushing and compaction of one story, minus the loss of gravity potential 
during the crushing of that story) exceeds the kinetic energy impacted to that story. Regardless of 
the load capacity of the columns, there is no way to deny the inevitability of progressive collapse 
driven by gravity alone if this criterion is satisfied (for the World Trade Center it is, with an order- 
of-magnitude margin). The parameters are the compaction ratio of a crushed story, the fracture of 
mass ejected outside the tower perimeter, and the energy dissipation per unit height. The last is the 
most important, yet the hardest to predict theoretically. Using inverse analysis, one could identify 
these parameters from a precise record of the motion of floors of a collapsing building. Due to a 
shroud of dust and smoke, the videos of WTC are useless here. It is proposed to obtain such records 
by monitoring the precise time history of displacements in different modes of building demolitions. 
The monitoring could be accomplished by real-time telemetry from sacrificial accelerometers, or 
by high-speed optical camera. The resulting information on energy absorption capability would be 
valuable for the rating of various structural systems and for inferring their collapse mode under 
extreme fire, internal explosion, external blast, impact or other kinds of terrorist attack, as well as 
earthquake and foundation movements. 

Introduction 

The destruction of the World Trade Center (WTC) on 9/11/01 was not only the biggest mass 
murder in the U.S. history but also a big surprise for the structural engineering profession, 
perhaps the biggest since the collapse of Tacoma Bridge in 1940. No experienced structural 
engineer watching the attack expected the WTC towers to collapse. No skyscraper has ever 
before collapsed due to fire. The fact that the WTC towers did, beckons deep examination. 

In this paper, 3 attention will be focused on the progressive collapse, triggered in WTC by 
fire and previously experienced in many tall buildings as a result of earthquake or explosions 
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(including terrorist attack). A simplified one-dimensional analytical solution of the collapse 
front propagation will be presented. It will be shown how this solution can be used to determine 
the energy absorption capability of individual stories if the motion history is precisely recorded. 
Because of the shroud of dust and smoke, these histories cannot be identified from the videos 
of the collapsing WTC towers, and so nothing can be learned in this regard from that collapse. 
However, monitoring of tall building demolitions, which represent a kind of progressive collapse, 
could provide such histories, and using different demolition modes can provide further valuable 
data. Development of a simple theory amenable to inverse analysis of these histories is the key. 
It would permit extracting valuable information on the energy absorption capability of various 
types of structural systems, and is, therefore, the main objective of this paper. 

Many disasters other than WTC attest to the danger of progressive collapse — e.g., the 
collapse of Ronan Point apartments in U.K., 1968 (Levy and Salvadori 1992), where kitchen 
gas explosion on the 18th floor sent a 25-story stack of rooms to the ground; the bombing 
of Murrah Federal Building in Oklahoma City, 1995, where the air blast pressure sufficed to 
take out only a few lower floors, while the upper floors failed by progressive collapse; the 2000 
Commonwealth Ave. tower in Boston, 1971, triggered by punching of insufficiently hardened 
slab; the New World Hotel in Singapore; many buildings in Armenian earthquake, Turkish, 
Mexico City and other earthquake; etc. A number of ancient towers failed in this way, too — 
e.g., the Civic Center of Pavia in 1989 (Binda et al. 1992); the cathedral in Goch, Germany; 
the Campanile in Venice in 1902; etc. (Heinle and Leonhardt 1989), where the trigger was 
centuries-long stress redistribution due to drying shrinkage and creep (Ferretti and Bazant 
2006). 

Review of Causes of WTC Collapse 

Although the structural damage inflicted by aircraft was severe, it was only local. Without 
stripping of much of the steel insulation during impact, the subsequent fire could not have led 
to overall collapse (Bazant and Zhou 2002, NIST 2005). As generally accepted in structural 
engineering and structural mechanics community (though not among some laymen seeking to 
unveil a conspiracy), the failure scenario was as follows: 

1. About 60% of the 60 columns of the impacted face of framed-tube (and about 13% 
of the total of 287 columns) were severed, and many more were significantly deflected. 
This caused stress redistribution, which significantly increased the load of some columns, 
attaining or nearing the load capacity for some of them. 

2. Significant amount of fire insulation was stripped during aircraft impact by flying debris 
(without that, the towers would likely have survived). In consequence, many structural 
steel members heated up to 600° C (NIST 2005) [the structural steel used loses about 
20% of its yield strength already at 300°C, and about 85% at 600°C, NIST 2005; and 
exhibits significant visco-plasticity, or creep, above 450° (e.g. Cottrell 1964, p. 299), 
especially in the columns overloaded by load redistribution; the press reports right after 
9/11, indicating temperature in excess of 800°C, turned out to be groundless, but Bazant 
and and Zhou's analysis did not depend on that]. 

3. Differential thermal expansion, combined with heat-induced viscoplastic deformation, 
caused the floor trusses to sag. The sagging trusses pulled many perimeter columns 
inward (by about 1 m, NIST 2005). The bowing of these columns served as a huge 
imperfection inducing multi-story out-of-plane buckling of framed tube wall. The lateral 



deflections of some columns due to aircraft impact and differential thermal expansion also 
decreased buckling strength. 

4. The combination of six effects — a) overload of some columns due to initial stress redistrib- 
ution, b) great lowering of yield limit and creep, c)lateral deflections of many columns due 
to thermal strains and sagging floor trusses, d) weakened lateral support due to reduced 
in-plane stiffness of sagging floors, e) multi-story buckling of some columns (for which the 
critical load is an order of magnitude less than it is for one-story buckling), and f) local 
plastic buckling of heated column webs — finally led to buckling of columns (Fig. lb). 
As a result, the upper part of tower fell, with little resistance, through at least one floor 
height, impacting the lower part of tower. This triggered progressive collapse because the 
kinetic energy of the falling upper part exceeded (by an order of magnitude) the energy 
that could be absorbed by limited plastic deformations and fracturing in the lower part 
of tower. 

In broad terms, this scenario was proposed by Bazant (2001), and Bazant and Zhou (2002) on 
the basis of simplified analysis relying on energy considerations. Up to the moment of collapse 
trigger, the foregoing scenario was identified by meticulous, exhaustive and very realistic com- 
puter simulations of unprecedented detail, conducted S. Shyam Sunder 's team at the National 
Institute of Standards and Technology (NIST). The progressive collapse was not simulated at 
NIST because its inevitability, once triggered by column buckling, had already been proven by 
Bazant and Zhou's (2002) comparison of kinetic energy to energy absorption capability. 

Before disappearing from view, the upper part of the South tower was seen to tilt signif- 
icantly (and of the North tower mildly). Some wondered why the tilting (Fig. Id) did not 
continue, so that the upper part would pivot about its base like a falling tree (see Fig. 4 of 
Bazant and Zhou 2002). However, such toppling to the side was impossible because the hori- 
zontal reaction to the rate of angular momentum of the upper part would have exceeded the 
elasto-plastic shear resistance of the story at least 10. 3x (Bazant and Zhou 2002). 

The kinetic energy of the top part of tower impacting the floor below was found to be 
about 8.4 x larger than the plastic energy absorption capability of the underlying story, and 
considerably higher than that if fracturing were taken into account (Bazant and Zhou 2002). 
This fact, along with the fact that, during the progressive collapse of underlying floors (Figs. Id 
and 2) the kinetic energy rapidly increases (roughly in proportion to the square of the number 
of stories traversed), sufficed to Bazant and Zhou (2002) to conclude that the tower was doomed 
once the top part of tower has dropped through the height of one story (or even 0.5 m). It 
was also observed that this conclusion made any calculations of the dynamics of progressive 
collapse after the first single-story drop of upper part superfluous. The relative smallness of 
energy absorption capability compared to the kinetic energy also sufficed to explain, without 
any further calculations, why the collapse duration could not have been much longer (say, twice 
as long or more) than the duration of a free fall from the tower top. 

Therefore, no further analysis has been necessary to prove that the WTC towers had to fall 
the way they did, due to gravity alone. However, a theory describing the progressive collapse 
dynamics beyond the initial trigger, with WTC as a paradigm, could nevertheless be very useful 
for other purposes, especially for learning from demolitions. Formulating this theory is the main 
objective of what follows. 



Motion of Crushing Columns of One Story and Energy Dissipation 

When the upper floor crashes into the lower one, with a layer of rubble between them, the 
initial height h of the story is reduced to \h, with A denoting the compaction ratio (in finite- 
strain theory, A is called the stretch). After that, the load can increase without bounds. In a 
one-dimensional model pursued here, one may use the estimate: 

A = (1 " KoutW/Vo (1) 

where Vo = initial volume of the tower, V\ ~ volume of the rubble on the ground into which the 
whole tower mass has been compacted, and K out = correction representing mainly the fraction 
of the rubble that has been ejected during collapse outside the perimeter of the tower and 
thus does not resist compaction. The rubble that has not been ejected during collapse but was 
pushed outside the tower perimeter only after landing on the heap on the ground should not be 
counted in K out . The volume of the rubble found outside the footprint of the tower, which can 
be measured by surveying the rubble heap on the ground after the collapse, is an upper bound 
on Vi, but probably much too high a bound for serving as an estimate. 

Let u denote the vertical displacement of the top floor relative to the floor below (Fig. 3, 4), 
and F(u) the corresponding vertical load that all the columns of the floor transmit. To analyze 
progressive collapse, the complete load-displacement diagram F(u) must be known (Fig. 3, 4 
top left). It begins by elastic shortening and, after the peak load, F steeply declines with u due 
to plastic buckling, combined with fracturing (for columns heated above approximately 450°C, 
the buckling is viscoplastic). For single column buckling, the inelastic deformation localizes 
into three plastic (or softening) hinges (see Figs. 2b, c and 5b in Bazant and Zhou 2002). For 
multi-story buckling, the load-deflection diagram has a similar shape but the ordinates are 
reduced roughly by an order of magnitude; in that case, the framed tube wall is likely to buckle 
as a plate, which requires four hinges to form on some columns lines and three on others (see 
Fig. 2c of Bazant and Zhou). Such a buckling mode is suggested by photographs of flying large 
fragments of the framed-tube wall, which show rows of what looks like broken-off plastic hinges. 

The mass of columns is assumed to be lumped, half and half, into the mass of the upper and 
lower floors, while the columns are massless. This avoids calculating the elastic waves propa- 
gating along the collapsing columns between the floors (which cannot interfere with buckling 
because their travel times are orders-of-magnitude shorter). 

Deceleration and Acceleration During the Crushing of One Story. The two in- 
tersections of the horizontal line F = gm(z) with the curve F(u) seen in Fig. 3 and 4a (top) 
are equilibrium states (there is also a third equilibrium state at intersection with the vertical 
line of rehardening upon contact). But any other state on this curve is a transient dynamic 
state, in which the difference from the line F = gm(z) represents the inertia force that must 
be generated by acceleration or deceleration of the block of tower mass m(z) above level z of 
the top floor of the story. 

Before impact by upper part, the columns are in equilibrium, i.e. F(u Q ) = gm(z), where 
Wo = initial elastic shortening of columns under weight gm(z) (about 0.0003/i). At impact, 
the initial condition for subsequent motion is velocity vq = u(uq) ~ Vi = velocity of the 
impacting block of upper part of tower. Precisely, from balance of linear momentum upon 
impact, v = m(z)/[m(z) + tuf], but this is only slightly less than Vi because rap <C m(z) (mj? 
is the mass of the impacted upper floor). 

When F(u) ^ gm(z), the difference F(u) — gm(z) causes deceleration of mass m z if positive 
(AFd in Fig. 3) and acceleration if negative (AF a in Fig. 3). The equation of motion of 
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mass m(z) during the crushing of one story (or one group of stories, in the case of multi-story 
buckling) reads as follows: 

u = g - F(u)/m(z) (2) 

where z = constant = coordinate of the top floor of the story, and superior dots denote deriv- 
atives with respect to time t. So, after impact, the column resistance causes mass m(z) to 
decelerate, but only until point u c at which the load-deflection diagram intersects the line 
F = gm(z) (Figs. 3, 4a). After that, mass m(z) accelerates until the end of column crushing. 

If the complete function F(u) is known, then calculation of the motion of the upper part 
of tower from (2) is easy (to calculate this function precisely is a formidable problem, but an 
upper bound curve is easy to figure out from plastic hinges, Bazant and Zhou 2002). Examples 
of accurately computed velocity v = u from Eq. (2) for various load-displacement diagrams 
graphically defined in the top row of Fig. 4a are shown in rows 2 and 3 of Fig. 4a,b,c. 

Energy Criterion of Progressive Collapse Trigger. The energy loss of the columns up 
to displacement u is 



$( u ) = / [F{v!) - gm{z)\ du = W{u) - gm(z)u (3) 

JUq 

W{u) = [ U F{u') du' (4) 
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where z = constant = column top coordinate, W(u) = energy dissipated by the columns = area 
under the load-displacement diagram (Fig. 3) and —gm(z)u = gravitational potential change 
causing an increment of kinetic energy of mass m(z). Note that, since the possibility of unload- 
ing (Fig. 4c top) can be dismissed, W(u) is path-independent and thus can be regarded, from 
the thermodynamic viewpoint, as the internal energy, or free energy, for very fast (adiabatic), 
or very slow (isothermal) deformations, and thus $(w) represents the potential energy loss. If 
F(u) < gm(z) for all u, $(«) continuously decreases. If not, then $(-u) first increases and then 
decreases during the collapse of each story. Clearly, collapse will get arrested if and only if the 
kinetic energy does not suffice for reaching the interval of accelerated motion, i.e., the interval 
of decreasing $(w), i.e., (Fig. 4 right column). So, the crushing of columns within one story 
will get arrested before completion (Fig. 4c) if and only if 

K, < W c (5) 

where W c = $(u c ) = W{u c ) — gm{z)u c = net energy loss up to u c during the crushing of one 
story, and K, = kinetic energy of the impacting mass m(z). This is the criterion of preventing 
progressive collapse to begin (Fig. 4c). Its violation triggers progressive collapse. 

Graphically, this criterion means that K, must be smaller than the area under the load- 
deflection diagram lying above the horizontal line F = gm(z) (Fig. 3 and 4 right column). If 
this condition is violated, the next story will again suffer an impact and the collapse process 
will be repeated. 

The next story will be impacted with higher kinetic energy if and only if 

W g > W p (6) 

where W g = gm(z)uf = loss of gravity when the upper part of tower is moved down by distance 
Uf, Uf = (1 — X)h = final displacement at full compaction, and W p = W{uf) = J f F(u)du 
= area under the complete load-displacement curve F(u) (Fig. 3). This is the criterion of 
accelerated collapse. 



For the case of WTC towers, it was estimated by Bazant and Zhou (2002) that /C ~ 
8AW P 3> W p for the story where progressive collapse initiated. W g was greater than W p by 
an order of magnitude, which guaranteed acceleration of collapse from one story to the next. 

Some investigators have been under the mistaken impression that collapse cannot occur 
if (because of safety factors used in design) the weight rag of the upper part is less than 
the load capacity F of the floor. This led them to postulate various strange ideas (such as 
"fracture wave" and planted explosives). However, the criterion in (5) makes it clear that this 
impression is erroneous. If this criterion is violated, there is (regardless of Fq) no way to deny 
the inevitability of progressive collapse driven only by gravity. 

Options for Transition to Global Continuum Model. One option would be finite 
element simulation based on the traditional homogenization of heterogenous microstructure 
of the tower, in which the load- displacement curve F(u) in Fig, 3 would be converted to an 
averaged stress-strain curve o~(e) by setting e = u/h and a = F/A) (A = cross-section area 
of the tower). However, the stress-strain relation delivered by this standard homogenization 
approach would exhibit strain softening. This would lead to an ill-posed dynamic problem, 
whose mathematical solution exists but is physically wrong (Bazant and Belytchko 1985, Bazant 
and Cedolin 2003, sec. 13.1). To obtain a well-posed formulation, it would be necessary to 
regularize the initial-boundary value problem by introducing a nonlocal formulation (Bazant 
and Jirasek 2004, Bazant and Cedolin 2003, chapter 13) with a characteristic length equal to 
the story height h (such regularization was forgotten in the "fracture wave" theory, proposed 
as an alternative for modeling WTC collapse]. But the nonlocal approach is relatively complex 
to program and gradual strain softening need not be modelled because only the total energy 
release per story is important (as demonstrated in rows 2 and 3 of Fig. 4 by equivalence of 
velocity diagrams). 

Because of the clearly defined periodic microstructure of stories in the tower, there is an- 
other, more effective, option: non-softening energetically equivalent characterization of discrete 
elements — the individual failing stories. This option is pursued next. It corresponds to non- 
standard homogenization, in which in which the aim is not homogenized stiffness but homoge- 
nized energy dissipation (this approach is analogous to the energy equivalent transition in van 
der Waals theory of gas-liquid phase changes, and the energy equivalence is also analogous to 
the crack band model for softening distributed damage, Bazant and Cedolin 2003, Bazant and 
Jirasek, 2004). 

Energetically Equivalent Mean Crushing Force. For the purpose of continuum smear- 
ing of a tower with many stories, the actual load-displacement diagram F(z) (curve OABC in 
Fig. 2a) can be replaced by a simple energetically equivalent diagram, represented by the hor- 
izontal line F = F c . Here F c is the mean crushing force (or resistance) at level z, such that 
the dissipated energy per story, represented by the rectangular area under the horizontal line 
F = F c , be equal to the total area W p under the actual load-displacement curve OABC; i.e., 

F c = ™> = if F(u) du (7) 



The energy-equivalent replacement avoids unstable snap-through (Bazant and Cedolin 2003) 
(and is analogous to what is in physics of phase transitions called the Maxwell line). Although 
the dynamic u(t) history for the replacement F c is not the same as for the actual F(u), the final 
values of displacement u and velocity u at the end of crushing of a story are exactly the same, 
as shown in the exactly calculated diagrams in rows 2 and 3 of Fig. 4. So the replacement has 
no effect on the overall change of velocity v of the collapsing story from the beginning to the 
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end of column crushing (Fig. 4), i.e, from u = to u = Uf (as long as F c is not large enough 
to arrest the downward motion). F c may also be regarded as the mean energy dissipated per 
unit height of the tower, which has the physical dimension of force. 

Note that it would be slightly more accurate not to include the minuscule elastic strain- 
energy portion of W p in integral (7), i.e., replace the lower limit with Uq. But then, instead of 
constant F c , we would need to consider an elastic-perfectly plastic force-displacement relation. 
The stiff and short elastic stress rise from u = to u (Fig. 4) would produce elastic waves 
which would create an unnecessary noise in the response and complicate numerical simulation). 

One- Dimensional Continuum Model for Crushing Front Propagation 

Detailed finite element analysis simulating plasticity and break-ups of all column and beams, 
and the flight and collisions of broken pieces, would be extremely difficult, as well as unsuited 
for extracting the basic general trends. Thus it appears reasonable to make four simplifying 
hypotheses: (i) The only displacements are vertical and only the mean of vertical displacement 
over the whole floor needs to be considered, (ii) Energy is dissipated only at the crushing front 
(this implies that the blocks in Fig. 2 may be treated rigid, i.e., deformations of the blocks 
away from the crushing front may be neglected), (iii) The relation of resisting normal force F 
transmitted by columns of each floor to the relative displacement u between two adjacent floors 
obeys a known load- displacement diagram (Fig. 4, terminating with a specified compaction 
ratio A (which must be adjusted to take into account lateral shedding of a certain known 
fraction of rubble outside the tower perimeter), (iv) The stories are so numerous, and the 
collapse front traverses so many stories, that a continuum smearing (i.e., homogenization) gives 
a sufficiently accurate overall picture. 

The one-dimensionally idealized progress of collapse of a tall building (of initial height H) 
is idealized in Fig. 2, where (, rj = coordinates measured from the initial and current tower top, 
respectively; z(t),y(t) = coordinates ( and r\ of the crushing front at time t (( is the Lagrangian 
coordinate of material points in the sense of finite strain theory, while y is measured from the 
moving top of building). The initial location of the first floor crashing into the one below is at 
Q = z = z = y . The resisting force F and compaction ratio A are known functions of z. A 
and C label the lower and upper undisturbed parts of tower, and B the zone of crushed stories 
compacted from initial thickness So to the current thickness 

r z(t) 

s(t) = / A(C)dC (8) 

When \x = constant, s(t) = X(z(t) — Zq) where z(t) — Zq is the distance that the crushing front 
has traversed through the tower at time t. The velocity of the upper part of tower is 

v(t) = [1 - \(z)]z(t) (9) 

First it needs to be decided whether crushed zone B will propagate down or up through the 
tower. The equation of motion of zone B requires that 

F ± - F 2 = \s [fig - (fiv)-} (10) 

This expression must be positive if the zone B is falling slower than a free fall, which is reasonable 
to expect and is confirmed by the solution to be given. Therefore always F 2 < F\. So, neither 
upward, nor two-sided simultaneous, propagation of crushing front is possible. 
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Note that a front propagating intermittently up and down would nevertheless be found 
possible if F c (z) were considered to be a random (autocorrelated) field. In that case, short 
intervals At may exist in which the difference F c \ — F C 2 of random F c values at bottom and top 
of crushed block B would exceed the right-hand side of Eq. (10). During those short intervals, 
crush-up would occur instead of crush-down, more frequently for a larger coefficient of variation. 
The greater the value of s , the larger the right-hand side of (10), and thus the smaller the 
chance of crush-up. So, random crush-up intervals could exist only at the beginning of collapse, 
when So is still small enough. Stochastic analysis, however, is beyond the scope of this study. 

The phase of downward propagation of the front will be called the crush-down phase, or 
Phase I (Fig. 4b). After the lower crushing front hits the ground, the upper crushing front of 
the compacted zone can begin propagating into the falling upper part of tower (Fig. 4d). This 
will be called the crush-up phase, or Phase II (it could also be called the 'demolition phase', 
because demolitions of buildings are usually effected by explosive cutter charges placed at the 
bottom). 

Let /* = //(C) = initial mass density at coordinate C = continuously smeared mass of undis- 
turbed tower per unit height. The mass density of the compacted zone B is mu(z)/X(z) (> fi). 
However, a correction must be made for the fraction K out of the mass that is being lost at the 
crushing front, ejected into the air outside the perimeter of the tower. During crash-down, the 
ejected mass alters the inertia and weight of the moving compacted part B, which requires a 
correction to m(z), while during crash- up no correction is needed because part B is not moving. 
Accordingly, we adjust the definition of the inertial mass of the tower above level z in the 
crush-down phase as follows: 

For z > z : m(z) = m(z ) + / (1 - /w)MC)dC, rn(z ) = / MC)dC (H) 

Jz JO 

No adjustment is needed for the crush-up phase because block B of compacted rubble does not 
move with C but is is stationary. 

Differential Equations of Progressive Collapse or Demolition 

The differential equations for z(t) and y(t) can be obtained from dynamic free body diagrams 
(Fig. 4h). In the crush-down phase, the compacted zone B and the upper part A of tower move 
together as one rigid body accreting mass, with combined momentum (1 — X)m(z)z. The neg- 
ative of the derivative of this momentum is the upward inertia force. Additional vertical forces 
are weight m(z)g downward, and resistance F c (z) upward. The condition of dynamic equilib- 
rium according to d'Alembert principle yields the following differential equation for compaction 
front propagation in the crush-down phase I of progressive collapse: 

— < m(z)[l — X(z)] — > — m(z)g = —F c (z) (crush-down) (12) 

(JLL I vj.6 I 

For the special case of A = F c = K out = 0, Eq. (12) reduces to (zz)' = z (the numerical solution 
for this special case was presented by Kausel, 2001). 

The initial conditions for the crush-down phase I are z = z and i = 0. Downward propa- 
gation will start if only if 

m(z )g > F c (z ) (13) 

In the crush-up phase, the crushing front at rj = y is moving up with velocity X(y)y, and 
so the downward momentum of part C is m(y)[l — X(y)]y. Downward acceleration of part C is 



opposed by upward inertia force 

F^ = -{rn(y)[l-\(y)]yy (14) 

By contrast to the crush-down phase, the compacted zone B with accreting mass is not moving 
with part C but is now stationary (Fig. 4d), and this makes a difference. During every time 
increment d£, the momentum 

dp = [fi(y)(ydt)][l - \(y)]y (15) 

of the infinitesimal slice dy = ydt at the crushing front gets reduced to (y < 0). So, the 
stationary part B is subjected to downward inertia force (Fig. 4g): 

F l B = dp/dt = f i(y)[l-\(y)]y 2 (16) 

(this is a similar phenomenon as, in the kinetic theory of gases, the pressure of gas molecules 
hitting a wall). As a reaction, the same force acts upward from part B onto part C. Adding 
also the force of gravity (and noting that y < 0,y < 0), the dynamic equilibrium of part C as 
a free body requires that F t B — Ff — m(y)g + F c = 0. This yields the following differential 
equation for compaction front propagation in the crush-up phase of progressive collapse: 




H-A< 



9 \ = F c (y) (crush-up) (17) 



For the special case of A = F c = and constant fi (for which m = jjy), Eq. (17) reduces to 
y = —g, which is the equation of free fall of a fixed mass. 

For the special case when only A is constant while F c (y) and ^(y) vary, Eq. (17) reduces to 

y = -g(y), ttv) = b - F c ( y )/m(y)]/(i - A) (is) 

This is equivalent to a fall under variable gravity acceleration g(y). Obviously, the collapse 
will accelerate (for A ^ 0) only as long as g > 0, i.e., if condition (13) is satisfied. Since 
\im y ^om(y) = 0, this condition will always become violated before collapse terminates (unless 
F c = 0), and so the collapse must decelerate at the end. 

For F c > 0, the tower can in fact never collapse totally, i.e, y = cannot be attained. To 
prove it, consider the opposite, i.e. y — > 0. Then y = C/y where C = F c /ji(l — A) = constant 
> 0; hence (y 2 )' = lyy = 2Cy/y, the integration of which gives y 2 = 2Cln(y/Ci) where C\ is a 
constant. The last equation cannot be satisfied for y — ► because the left-hand side > while 
the right-hand side — ► — oo; Q.E.D. 

As the rubble height approaches its final value, i.e. for lim^ = 2//(> 0)> the values of 
m,X,F c are nearly constant, and so y = (F c /m — g)/(l — A) = Cq = constant [> 0, which is 
again condition (13]. Hence, y = Co, which gives y{t) — yf = Co(t — t/) 2 . So, if F c > 0, the 
collapse history y(t) will terminate asymptotically as a parabola at some finite height y\ and 
finite time tf. 

For a more detailed simulation of collapse, it would be possible to use for each floor Eq. 
(2) for motion within each floor, or introduce into Eq. (12) and (17) a function F c (z) varying 
within each story height as shown by the actual response curves in Figs. 4 and 5. This would 
give a fluctuating response with oscillations superposed on the same mean trend of z(t) or y(t) 
as that for smooth F c (z). Extremely small time steps would be needed in this case. 



The fact that F c is smaller in the heated story than in the cold stories may be taken into 
account by reducing F c (z) within a certain interval z G (zq,Z\). 

The initial conditions for the crush-up phase II are y = y = z and a velocity y equal to the 
terminal velocity of the crush-down phase. For a demolition, triggered at the base of building, 
the initial conditions are y = y and y = 0, while F c = for the y value corresponding to the 
ground story height. 

For trigger by an explosion or vertical impact, the present formulation might be used with an 
initial condition consisting of a certain finite initial velocity Vq caused by the explosion. In that 
case, /C in collapse trigger criterion (5)may be replaced by energy imparted by the explosion. 

Dimensionless Formulation 

To convert the formulation to a dimensionless form, note that the solution can be considered to 
be a function of 2 coordinates, t and z (or y), and 6 independent parameters, H, Zo,g, F c , fi(z), \(z), 
and involves involves 3 independent dimensions, the mass, length and time. According to the 
Vashy-Buckingham theorem, the solution must depend on only 7 + 2 — 3 = 6 dimensionless 
independent parameters, of which 2 are the dimensionless time and spatial coordinate. They 
may be chosen as follows: 



r = t y/g/H, Z = z/H or Y = y/H, Z = z /H = y /H 

F C (Z) = F c (z)/Mg, m(Z) = m(z)/M, A = \{Z) (19) 

where M = m(H) = total mass of the tower. After transformation to these variables, the 
differential equations of the problem, Eqs. (12) and (17), take the following dimensionless 
forms: 

r\ ( dZ) 

[1 - \{Z)\ fh(Z) — \ - fh(Z) = -F C {Z) (crush-down) (20) 



dr I dr 



fh(Y) 




F C (Y) (crush-up) (21) 



The dimensionless form of the initial conditions is obvious. 

In the special case of constant [i and A, we have m(Z) = Z, m(Y) = Y , and the foregoing 
dimensionless differential equations take the form: 

(1- \)(ZZ + Z 2 )-Z = -F C (Z) (crush-down) (22) 

(1 - \)YY + Y = F C {Y) (crush-up) (23) 

Numerical Solution and Parametric Study 

Eq. (12) may be converted to a system of 2 first-order differential equations of the form z = x 
and x = F(x, z), with prescribed values of z and x as the initial conditions. This system can be 
easily solved by some efficient standard numerical algorithm, such as the Runge-Kutta method. 
The same is true for Eq. (17). 

The diagrams in Fig. 6 present the collapse histories computed for the approximate parame- 
ters of WTC (heavy solid curves) and for modified values of these parameters. For comparison, 
the curve of free fall from the tower top is shown in each diagram as the leftmost curve. The 
transition from the crush-down phase I to the crush-up phase II is marked in each diagram 
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(except one) by a horizontal line. The parameter values used for calculation are listed in each 
diagram. There were chosen as the typical values for the WTC and their variations. Wf denotes 
the mean of a linearly varying crushing energy Wf. Since the first story to collapse was heated, 
the value of F c within the interval of z corresponding to the height of that story was reduced to 
one half. Fig. 7 shows separately the histories of tower top coordinate for the crush-up phase 
alone, which is the case of demolition. Four characteristics of the plots of numerical results in 
Fig. 6 and 7 should be noticed: 

• Varying the building characteristics, particularly the crushing energy Wf per story, makes 
a large enough difference in response to be easily detectable by the monitoring of collapse. 

• The effect of crushing energy Wf on the rate of progressive collapse is much higher than 
the effect of compaction ratio A or specific mass \x. This means that these two parameters 
need not be estimated very accurately in advance of inverse analysis. 

• For the structural system such as WTC, the energy required to arrest the collapse after 
a drop of only one or several stories (Fig. 6e) would have to be an order of magnitude 
higher than the energy dissipation capacity of the structural system used in WTC . 

• For the typical WTC characteristics, the collapse takes about 10.8 s (Fig. 6 top left), 
which is not much longer (precisely only 17% longer) than the duration of free fall in 
vacuum from tower top to the ground, which is 9.21 s (the duration of 10.8 s is within 
the range of Bazant and Zhou's, 2002, crude estimate). For all of the wide range of 
parameter values considered in Fig. 6, the collapse takes less than about the double of 
free fall duration. 

The last two points confirm Bazant and Zhou's (2002) observations about collapse duration 
made on the basis of initial kinetic energy and without any calculation of collapse history. 

What Can We Learn? — Proposal for Monitoring Demolitions 

We have seen that the main unknown in predicting cohesive collapse is the mean energy dis- 
sipation Wf per story. The variable fi(z) is known from the design, and the contraction ratio 
X(z) can be reasonably estimated from Eq. (1) based on observing the rubble heap after col- 
lapse. But a theoretical or computational prediction of F c is extremely difficult and fraught 
with uncertainty. Precise collapse observations are required. 

Eqs. (12) and (17) show that F c (z) can be evaluated from the monitoring of motion history 
z(t) and y(t), provided that ji(z) and X(z) are known. Such information can, in theory, be 
extracted from a high-speed camera record of the collapse. Approximate information could 
have been extracted from a regular video of collapse if the moving part of the towers of WTC 
were not shrouded in a cloud of dust and smoke. Thus, despite thousands of photographic 
records, nothing can be learned from WTC. 

However, valuable information on the energy dissipation capability of various types of struc- 
tural systems could be extracted by monitoring demolitions. During the initial period of de- 
molition, the precise history of motion of building top could be determined from a video of 
collapse. After the whole building disappears in dust cloud, various remote sensing techniques 
could be used. For example, one could follow through the dust cloud the motion of sacrifi- 
cial radio transmitters. Or one could install sacrificial accelerometers monitored by real-time 
telemetry. From the acceleration record, the y(t)-history could be integrated. 
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Therefore, monitoring of demolitions is proposed as a means of learning about the energy 
absorption capability of various structural systems. 

Usefulness of Varying Demolition Mode. Ronan Point apartments, the Oklahoma City 
bombing, etc., demonstrate that only a vertical slice of building may undergo progressive col- 
lapse, while the remainder of building stands. Such a collapse is truly a three-dimensional prob- 
lem, much harder to analyze, but some cases might allow adapting the present one-dimensional 
model as an approximation. For example, in Ronan Point apartments, energy was dissipated 
not only by vertical crushing of stories, but also by shearing successive floor slabs from their 
attachments to columns on the side of the collapsing stack of rooms. The present model seems 
usable if the energy dissipated by shearing is added to the crushing energy F c , and if the rota- 
tional kinetic energy of floor slabs whose fall is hindered on one side by column attachments is 
taken into account. Such a generalization of the present model could be calibrated by compar- 
ing data from two different demolition modes: 1) the usual mode, in which the building is made 
to collapse symmetrically, and 2) another mode in which first only a vertical slice of building 
(e.g., one stack of rooms) is made to collapse by asymmetrically placed cutter charges. Many 
variants of this kind may be worth studying. 

Complex Three-Dimensional Situations. Situations such as stepped tall buildings call 
for three-dimensional analysis. Large-scale finite-strain computer simulation tracking the con- 
tacts of all the pieces of crushing floors and columns could in principle do the job but would 
be extraordinarily tedious to program and computationally demanding. The present analysis 
would be useful for calibrating such a computer program. 

Massive Structures. Progressive collapse is not out of question even for the massive load- 
bearing concrete cores of the tallest recent skyscrapers, as well as for tall bridge piers and 
tall towers of suspension or cable-stayed bridges (that such a collapse mode is a possibility is 
documented, e.g., by the collapses of Campanile in Venice and Civic Center tower in Pavia). 
Although progressive collapse of the modern massive piers and towers would be much harder 
to initiate, a terrorist attack of sufficient magnitude might not be inconceivable. Once a local 
damage causes a sufficient downward movement of the superior part of structure, it cannot be 
stopped. The question is, for instance, whether it might be within the means of a terrorist to 
cause, e.g., formation and slipping of an inclined crack band typical of compression fracture of 
concrete. In this regard, note that the size effect in compression fracture (Cusatis and Bazant 
2006) would assist a terrorist. 

Alternative Formulations, Extensions, Ramifications 

Alternative Derivation. A more elementary way to derive the differential equation for the 
crush- up phase is to calculate first the normal force N(rj) (positive if tensile) in cross section 
of any coordinate r\ G (0,y) (Fig. 4h). The downward velocity of block C is v = [1 — X(y)]y, 
and its acceleration is opposed by inertia force [1 — \{y)]ym{rj). The downward gravity force 
on this block is gm{rj). From dynamic equilibrium, the normal force N{rj) (positive if tensile), 
acting at the lower face rj of this block, is: 

N(rf) = -[I - \{y)\ y m{rj) + g m{r}) (24) 

For the crushing front, rj = y, this must be equal to the crushing force, i.e., N(y) = —F c (y). 
This immediately verifies Eq. 17. 

For the crush-down phase, the same expression holds for the cross section force N(Q. How- 
ever, in the dynamic equilibrium condition of block C, one must add upward inertia force /d(z)z 2 
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needed to accelerate from to i the mass that is accreting to block C per unit time. This then 
verifies Eq. 12. 

Potential and Kinetic Energies. An energy based formulation is useful for various ap- 
proximations, numerical algorithms and bounds. It is slightly complicated by the accretion of 
mass to the moving block and the dissipation of energy by crushing force F c . 

Consider first the crush-down phase. Since unloading of columns does not occur, a potential 
II can be defined as the gravitational potential minus the work of F c . Its rate is: 

^- = {F c [z(t)\ ~ 9m[z(t)]}v(t) (25) 

Due to accretion of mass to the moving block, its kinetic energy m(z)v 2 /2 in increased by 
the kinetic energy due to accelerating every infinitesimal slice dz = zdt of mass m'(z)(zdt) to 
velocity v. This means that kinetic energy increment \[m! {z){zdt)\v 2 is added during every 
time increment dt. So, the rate of added kinetic energy is ^m'(z)zv 2 , and the overall rate of 
change of kinetic energy /C is 

d/C(t) d fl [n]2M | 1 ,, 2 , ,dz(t) 



dt dt 12 



m[z(t)]v 2 (t)} + -m'(z)v 2 (t)^- (26) 



where m'(z) = dm(z)/dz, which differs from ji(z) only at locations where n out has an effect. 
Conservation of energy requires the sum of the last two energy rates to vanish. This condition 
yields: 

m(z)vv + \rri \z)(z)v 2 'z + \m '{z){z)zv 2 + [F c {z) — gm(z)]v = (27) 

Dividing this equation by mass velocity v and setting v = (1 — \)z, we find that Eq. (17) 
ensues. This verifies correctness of the foregoing energy expressions for the crush-down phase. 
For the crush-up phase, the rate of energy potential is 

^ = {9m[y(t)] - F c [y(t)]}v(t) (28) 

In formulating the kinetic energy, there is a difference from crush-down: The mass of each 
infinitesimal slice dy = ydt is during dt decelerated from velocity v to 0, removed from the 
moving block C, and added to the stationary block B. By analogous reasoning, one gets for the 
kinetic energy rate the expression: 

dK( " - d <\ mly{t) y (t) \-' m y {t m (29) 



dt dt 12 L * v /J v ') T vay /J v ' dt 

where ji(y) = m'(y). Energy conservation dictates that the sum of the last two energy rate 
expressions must vanish, and so 

m(y)vv + \n(y)v 2 y - \n(y)yv 2 + [gm(z) - F c (z)}v = (30) 

After division by v = (1 — X)y, Eq. (12) for the crush-up phase is recovered. This agreement 
verifies the correctness of the foregoing energy rate expressions. 

Lagrange equations of motion or Hamilton's principle (Flugge 1962) are often the best way to 
analyze complex dynamic systems. So why hasn't this approach been followed? — Because these 
equations are not valid for systems with variable mass. Although various special extensions to 
such systems have been formulated, they are complicated and depend on the particular type of 
system (e.g., Pesce 2003). 
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Solution by Quadratures for Constant A and fi, and K out = 0. In this case, which 
may serve as a test case for finite element program, Eq. (12) for the crush-down phase takes 
the form: 

ff + f 2 -Qf = -P or (ffY = Qf-P (31) 

Here Q = 1/(1 - A), P(t) = F c /fi(l - X)gH, F c = F c [z(t)), f = f(t) = z{t)/H\ and the superior 



dots now denote derivatives with respect to dimensionless time r = tJg/H. Let <p = f 2 /2. 
Then ip = ff and 



(p = Q^/2cp-P (32) 

(V? 2 )'=2^ = 2(Q v / 2^-P) V 9 (33) 

d(<p 2 ) = [ 2(Qjty - P)d<p (34) 



p> = (±QV2^ 3/2 -2P ¥ > + c) 1/2 (35) 



To 



r ' ( lQV2(p 3/2 - 2Pp) + C) d<p (36) 



><p(to) 



The second equation was obtained by multiplying the first by 2p>, and Eq. 35 was integrated 
by separation of variables; C and To are integration constants defined by the initial conditions. 
The last equation describes the collapse history parametrically; for any chosen (p, it yields the 



time as t = zJH/g or yJH/g where z or y = Hy/2ip. 

Eq. (12) for the crush-up phase with constant fi and A takes the form: 

ff + Qf = P (37) 

Multiplying this equation by /// and noting that // = \{f 2 )' and /// = (In/)', one may get 
the solution as follows: 

(f 2 Y = 2(Pf/f-Qf) (38) 

f 2 = 2(Plnf-Qf) + C (39) 

d/ = [2(P\nf-Qf) + C] l ' 2 &T (40) 

/ /(r) [2(Pln/-Q/) + C]- 1 / 2 dr (41) 

•'/(to) 



To 



Effect of Elastic Waves. Note that Eqs. (12) and (17) do not represent elastic stiffness 
and thus do not model wave propagation. Elastic waves are ignored because the maximum 
shortening of story height for which the columns remain elastic is only about 0.4 mm. If elastic 
response were incorporated, a much faster elastic wave, with step wavefront equal to Fq, would 
be found to emanate from the crushing front when each floor is hit, propagate down the tower, 
reflect from the ground, etc. But the damage to the tower would be almost nil, the subsequent 
propagation of the crushing front would be almost unaffected by the precursor elastic waves, 
and the solution would get more complicated. 

Analogous problem — Crushing of Foam. One-dimensional crushing of rigid foam by 
impact can be solved from the present differential equation for the crush-down phase. To this 
end, the top part of tower needs to be replaced by a rigid impacting object of mass equivalent 
to m(z ), whose initial velocity is assigned as the initial condition at t = 0. Compared to inertia 
forces, normally one can neglect gravity (i.e., set g = 0). 
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Implications and Conclusions 

1. If the total (internal) energy loss during the crushing of one story (representing the energy 
dissipated by the complete crushing and compaction of one story, minus the loss of gravity 
potential during the crushing of that story), exceeds the kinetic energy impacted to that 
story, collapse will continue to the next story This is the criterion of progressive collapse 
trigger (Eq. 5). If it is satisfied, there is no way to deny the inevitability of progressive 
collapse driven by gravity alone (regardless of by how much the combined strength of 
columns of one floor may exceed the weight of the part of tower above that floor). 

2. One-dimensional continuum idealization of progressive collapse is amenable to a simple 
analytical solution which brings to light the salient properties of the collapse process. 
The key idea is not to use classical homogenization, leading to a softening stress-strain 
relation necessitating nonlocal finite element analysis, but to characterize the stories by 
an energetically equivalent snapthrough. 

3. Distinction must be made between crush-down and crush-up phases, for which the crush- 
ing front of a moving block with accreting mass propagates into the stationary stories 
below, or into the moving stories above, respectively. This leads to a second-order non- 
linear differential equation for propagation of the crushing front in the crush-down phase, 
or in the subsequent crush-up phase. 

4. The mode and duration of collapse of WTC towers are consistent with the derived model, 
but nothing more can be learned because the motion was obstructed from view by a cloud 
of dust and smoke. 

5. The present idealized model allows simple inverse analysis which can yield the crushing 
energy per story and other properties of the structure from a precisely recorded history of 
motion during collapse. From the crushing energy, one can infer the collapse mode, e.g., 
single-story or multi-story buckling of columns. 

6. It is proposed to monitor the precise time history of displacements in building demolitions— 
for example by radio telemetry from sacrificial accelerometers, or high-speed optical 
camera — and to engineer different modes of collapse to be monitored. This should pro- 
vide invaluable information on energy absorption capability of various structural systems, 
needed for assessing the effects of explosions, impacts, earthquake, and terrorist acts. 
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